clear
clc;
clf;
L1=4000;
L2=4000;
V1 = 82.5;
V2 = 100;
% daofun1 = @(t,y)fly_the_air(t,y,L1,L2,V1,V2,0.5);
t= 0:1:360;
[t,xxyy] = ode45( @(t,y)fly_the_air(t,y,L1,L2,V1,V2,0.2),t,[10000,45000,50000,50000,50000,20000]);

plot(xxyy(:,1),xxyy(:,2),'k');hold on;
plot(xxyy(:,3),xxyy(:,4),'b');hold on ;
plot(xxyy(:,5),xxyy(:,6),'r');hold on;
axis([0 50000 0 70000])